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Summary 

Arcing to negatively biased, exposed solar cell 
interconnects on solar arrays placed in plasma 
environments has been well established in laboratory tests 
and inferred from space data. Such arcing may cause 
damaging interference with the operation of electrical 
power systems in spacecraft planned to be driven with 
high-voltage solar arrays. A simple analytical model has 
been developed to estimate the effects of negative bias 
arcs on solar array power system performance. Solar cell 
characteristics, plasma interactions, and power system 
features are modeled approximately by a linear, lumped 
element transient circuit, and the time domain equations 
are solved. Exact numerical results for solar array 
common-mode and load voltage transients are calculated 
for typical conditions. Acceptable load transients are 
found for a range of arc current amplitudes and time 
constants. 


Introduction 

Laboratory and flight tests of small, externally biased 
solar cell arrays show that arcs occur on or near 
negatively biased solar cell interconnects exposed to 
plasmas typical of low Earth orbit (LEO) (refs. 1 to 4). 
The voltage threshold of these arcs appears to lie several 
hundred volts below the plasma ground potential. The 
large solar arrays desired for future space missions in 
LEO will have areas operating at voltages more negative 
than this threshold, if their terminal voltages are 
sufficiently high (ref. 5). This is especially true if the 
arrays are permitted to float electrically in the space 
plasma. Negative bias arcing is a threat to the operation 
of such arrays if it can be shown that unacceptable 
electrical transients in the power system due to arcing 
cannot be avoided. Conversely, arcing may be tolerated if 
transients and other arc effects are found to be or can be 
made minimal. 

In order to estimate the severity of transients produced 
by negative bias arcs, a circuit model of a simple solar 
array power system has been developed. It shows 
approximately how load transients due to arcs on an 
array in LEO depend on plasma parameters, stored 
charge, solar cell behavior, and arc characteristics. 


Suitable linearizations permit a completely analytical 
solution to be obtained in the time domain. 

The model is current driven. Current pulses are 
injected at the negative end of an illuminated and loaded 
solar array that is in prior electrical equilibrium with the 
plasma. A pulse causes a positive, common-mode voltage 
transient of the array circuit relative to the plasma and 
perturbs the solar cell currents and voltages so as to 
produce a transient across the load. The size of this 
transient depends in part on the operating points of the 
solar cells. In this report we assume operation on the 
high-voltage side of the maximum power point where the 
internal dynamic resistance of the cells is low. Operation 
at other points can be accommodated by the model but is 
expected to result in much larger voltage transients (see 
the Solar Cell and Solar Array Characteristics and the 
Results and Discussion sections.) 

Several features of the model establish it as a worst- 
case analysis. First, the arc injection point is placed at the 
negative end of the array where it is expected, 
theoretically, to cause larger transients than if placed at 
more positive points. Second, the spatial distribution of 
the plasma interaction w'ith the array is approximated by 
lumped circuit elements arranged to allow upper-bound 
estimates of load transients. Third, the load is assumed to 
be purely resistive. The filtering expected to be included 
in actual power system loads is withheld. These features 
permit a rather simple model to be made, one that is 
useful for probing the space of arc current parameters 
yielding acceptable load transients on a given array under 
various conditions. 

Solar Cell and Solar Array 
Characteristics 

A solar array consisting of M parallel chains of 
matched solar cells is assumed. Each chain contains N 
cells wired in series, and each chain is connected at its 
ends to other chains. Arcs injected at one end of the 
chains cause equal voltage distributions along them. 
Therefore, cross wiring between adjacent points of the 
chains is omitted. Furthermore, model linearity requires 
that the solar cells not be permitted to depart greatly 
from their quiescent operating points during a transient. 
Limiting the arc current amplitude achieves this and 
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makes modeling of reverse biasing, short- and open- 
circuiting of solar cells, and diode protection 
unnecessary. 

If the array is illuminated and loaded so that each cell 
operates at quiescent voltage and current (e 0 and / 0 ), then 
the terminal voltage F^is Ne (J , the output current I () is 
Mi ()y the output power is MNe 0 i 0 , and the load resistance 
R l is Nt 0 /Mi 0 . For small fluctuations about the 
operating point, the terminal voltage is given by 

V 7 — n-PaV-la) 0 ) 

where p () ~Np/M and p is the magnitude of the dynamic 
resistance dt/di at the operating point of a single cell 
having current /voltage characteristic <(/). In the model, 
transient currents are confined to those for which 
equation (1) is a good approximation. 

The characteristic for a silicon solar cell may be 
approximated by the function 

c-(/) = </(ln (/j>-/)) + c 

where, for a 2 cm by 2 cm cell at normal operating 
temperature under full solar illumination, the values 
a = 40 mV, h- 132 mA, and c = 355 mV apply (ref. 6). At 
full power, an array of such cells having N~ 1000 and 
M = 1 500 produces 82 kW at 450 V with R l = p 0 = 2.5 
ohms. At the high-voltage 80-percent power point, which 
is a typical operating point used for model calculations, 
tiie array produces 66 kW of power at 508 V with R t - 4 
ohms and = 0.6 ohm. Operation on the high current 
side of maximum power should be avoided in real power 
systems if large arcs are expected. (See the Results and 
Discussion section.) 

A silicon solar cell has an internal capacitance that 
appears in parallel with p. For 2 cm by 2 cm silicon cells, 
C'0.02 p\\ Thus, the impedance due to C becomes 
important, as compared with p, only for pulses lasting 
about a microsecond or less. Although such short pulses 
are seen infrequently in laboratory tests, pulses with very 
short rise times do occur. Therefore, C must be included 
to model their responses correctly. 

Plasma Interactions 

Plasma Current Collection 

T he plasma environment used in this report is that of 
the ionosphere a near 400-kni altitude (ref. 7). Plasma 
densities in this region rise to 10 s or l() (> cm \ and 
electron and ion mean thermal energies are a few tenths 
of an electron volt. The ions are () * , predominantly. 
Both the Debye length and the electron gymradius are 


greater than a centimeter. The electron plasma frequency 
exceeds 1 MHz. 

A solar array placed in such an environment draws 
current from the plasma to the exposed metallic 
interconnects. Inasmuch as the interconnects are small 
compared to the Debye length and electron gyroradius, 
one might expect orbit-limited Langmuir probe theory to 
give an adequate estimation of current collection from 
the plasma. For electrons this theory predicts currents to 
a sphere at voltage V to be proportional to the collecting 
area, the electron thermal current density, and the factor 
(1 +eV/kT c ). Ion collection in this theory is somewhat 
more involved (ref. 8). 

Experimentally the situation is not so simple. Whereas 
currents to negatively biased interconnects appear to vary 
roughly linearly with voltage, those to positively biased 
interconnects are highly nonlinear. They vary as 
(V- V () Y ? where /?= 1/2 at low voltage, p> 1 in a snap- 
over region (ref. 1) where secondary emission effects on 
the cover slides become important, and p- 1 at high 
voltage. However, since snap-over is known 
experimentally to require seconds to develop, it is 
unlikely to play a role in the plasma response to the 
micro- and millisecond arcs expected. Thus, the voltage 
dependence of electron collection may be much more 
slowly varying on these time scales than is indicated by 
the quasistatic measurements in which snap-over is seen. 
The exact dependence is not yet well understood 
experimentally. 

To make progress toward an analytical model, 
simplifying assumptions regarding plasma current 
collection must be made. The most helpful is that both 
electron and ion plasma currents are ohmic, i.e., 
proportional to voltage. Such currents are characterized 
by conductances that depend on the properties of the 
ambient plasma and on the collecting area per cell. This is 
the equivalent of orbit-limited, spherical Langmuir probe 
collection in a thermal plasma with eV>>kT. Since 
eV> > k 7 almost everywhere on a high voltage array in 
LEO, the ohmic assumption treats ion collection 
realistically but ignores electron snap-over. 

We will assume for species s (electrons or () 1 ions) that 

I s = tx s V (2) 

where 




The plasma is assumed to be Maxwellian; n s is the particle 
density; T s and m s are the temperature and particle mass; 
A is the collecting area per cell; and j s is the ambient 
thermal current density. Typical values of the parameters 
used in the model are n e -n { - 4xl0 5 cm^ 3 , kT e ~0.2 
eV - 2kT h and A = 1 cm 2 . Then, a e = 2A x 10 ~ 6 mho for 
the electrons, and ai = 2xl0~ 8 mho for the ions. 

In solar array orientations where ram ion currents to 
negative interconnects dominate the thermal ion currents, 
the previous formulation requires adjustment. It may be 
assumed that ion collection is then given by spherical 
Langmuir probe theory in a streaming plasma in the limit 
where eV>>E 0 (ion ram energy). Thus, 

I r = a r V (3) 

where 

a r = e 2 Aiij( — L- ) 1/2 
\ m iE 0 / 

The ram ion energy for 0+ at 400 km is about 5 eV. This 
value gives a r ~ 10 - 8 mho. Thus, for this case, ram ion 
collection is not radically different in magnitude than 
thermal ion collection, and numerical results for this case 
will not be given. 

Stored Charge 

The several parts of a solar cell in an array are 
eapacitively coupled to each other and to the surrounding 
plasma. Substrate to cover slide, cover slide to plasma 
sheath, and interconnect to sheath capacitances in 
series/parallel combination form an equivalent 
capacitance per cell (C p ). An estimate of this capacitance 
lies in the range 1 to 200 pF for standard 2 cm by 2 cm 
silicon cells. The low value is approximately the 
capacitance to space of a sphere having area equal to that 
of one interconnect. The high value is the capacitance of 
a 2 cm by 2 cm cover slide to its substrate. Which of these 
limits dominates C p depends on whether a cover slide can 
maintain itself near the plasma potential during an arc by 
charging from the plasma. If so, C p tends toward the 
higher limit. Because experimental values are not yet well 
known, the present calculations draw upon the whole 
range of C p . 

Floating Equilibrium 

On a self-biased solar array in electrical equilibrium 
with the space plasma, the electron versus ion current 
balance condition determines what fraction of the array 
lies above the plasma potential and what fraction below. 
Those areas above the plasma potential collect electrons 
and those below collect ions, almost exclusively. In the 
ohmic collection approximation, the conductances a e and 


otj determine how the array floats electrically. For 
a e = 2Ax 10-6 m ho and a,- = 2x 10” 8 mho an array will 
float with about 8 percent of its area positive (see 
appendix A). 

During an arc, while individual cells are being driven 
positive in voltage, it is expected that a cell’s conductance 
jumps from a/ to a e as the cell passes the plasma potential 
and back again on the negative swing. Unfortunately, 
such voltage-dependent conductance switches are beyond 
the scope of the present work. In this model, once 
floating equilibrium is established and all conductances 
determined prior to an arc, they are not permitted to 
change. It is possible, however, to partially explore the 
consequences of conductance switching and the resulting 
enhanced electron collection during an arc by using 
electron conductances somewhat greater than those 
expected physically. 

Equivalent Circuits 

Transient Response of a Chain of Cells 

Figure 1 shows a transient circuit for a single chain of 
solar cells. Z p is the impedance of the parallel 
combination of solar cell plasma capacitance C p and 
resistance 1/a, and p and C have their previous 
definitions. An arc current injected at K_ flows toward 
K + and the load while leaking out to the plasma sheath 
along the way. The load transient is developed as the sum 
of the voltage transients across the cells due to the current 
transients through them. 

Since current flowing in a chain toward V + is 
diminished by leakage to the sheath through the Z p s, cell 
voltage transients also diminish in this direction. Thus, if 
all the Z p s were placed in parallel at V + , an overestimate 
of the load transient would result because all the cells 
would see the maximum current. Similarly, placing all the 
Z p s at would result in no load transient. These two 
cases bound the actual, uniform distribution case. An 
appropriate mixture of these cases approximates the 
actual one. 

Using a mixture of the bounding cases results in a 
considerable simplification of the plasma interaction 
circuits. The Z p s in each chain are split into two groups, 


c c c c c 





1 Plasma sheath 

figure 1.— Cell chain transient circuit. 
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one placed at K_ and the other at V + , and the groups are 
combined in parallel at each terminal. A plasma split 
parameter 6 adjusts the relative numbers of Z p s at each 
terminal while keeping the whole constant (see appendix 
A). Ranging this parameter permits continuous variation 
from one bounding case to the other. Figure 2 shows the 
mixed-case circuit for a single chain of cells. 

Solar Array Equivalent Circuit 

Figure 3 shows the resultant solar array power system 
transient circuit used in the model. /? + ,/?_,(?+, and C_ 
are plasma interaction components whose values are 
determined as described in appendix A. C a = MC/N and 
p G and Ri have their previous definitions (see the Solar 
Cell and Solar Array Characteristics section). 

As shown in appendix B, the terminal transient 
voltages, v+ and v , satisfy a set of coupled, first-order 
differential equations with the arc pulse i s (t) as current 
driver. The equations can be solved for arbitrary i s (t) by 
the method of Laplace transforms. With 

‘s(o = h [ ex p ( -y ) ~ exp ( ~ir ) 


the time-domain solutions have the form 
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Figure 2.— Mixed-case single chain. 
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Figure 3. —Array transient circuit. 


v(t)=B l exp ^ ~y ~ ) +Bl exp {~f~) +Bi 


+ B 4 exp (X 2 /) 


( 4 ) 


where B x ... B 4 , X b and X 2 are functions of the circuit 
parameters and of / 0 , K , and T. With K> 1, the driving 
current pulses have approximately the same time 
dependence as do the arcs seen in laboratory tests. 

Although an alternate parameterization of the arcs 
using an R-L-C series circuit is possible, the l ()y T, and K 
parameterization is more directly comparable to 
experimental measurements. 


Results and Discussion 

Figures 4 to 8 show typical numerical results produced 
by the model. In the figures, V- is the voltage transient 
amplitude at the negative end of the array and u ( j is the 
load transient amplitude. The values N = 1000, M - 1500, 
K = 5, a e = 2.4x10 mho, and a/ = 2.10-« mho apply 
throughout. Q is the total charge released, i.e., the 
integral of i s (t). The arc current amplitude I 0 is limited to 
200 A to insure the validity of the solar cell characteristic 
linearization (see the Solar Cell and Solar Array 
Characteristics section). 

Figure 4 shows the first-quadrant solar cell 
characteristic curve used. Operating points for various 
fractions of maximum power are shown. The sharply 
increasing negative slope to the right of the maximum 
power operating point suggests that transient amplitudes 
per unit arc current should be larger in this region than 
they are at operating points to the left of the maximum 
power point where the slope and, hence, the dynamic 



Figure 4. — Solar cell characteristic. Open-circuit voltage, < * , 550 mV; 
short-circuit current, / Sl , 132 mA; voltage at maximum power, 

450 mV; current at maximum power, 122 in A; maximum 

power, P lmn , 55 111 W. 
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Figure 5. — Transient amplitudes as function of maximum power 
fraction. 


resistance is smaller. Figure 5, which shows transient 
voltage amplitudes as a function of operating power for 
typical conditions, bears this out. There is a factor of 
almost three between values of v j at the 90-percent points 
on either side of maximum power. This shows that 
transient amplitudes can be reduced significantly by 
operating on the high-voltage side of maximum power. 

Amplitudes shown in figures 6 to 8 have been 
calculated assuming operation at 80-percent power in the 
high-voltage side. Figure 4 shows that solar cell 
linearization is a good approximation at this operating 
point if cell current transients are kept below about 30 
mA; / o <200 A achieves this. 



Figure 6. — Transient amplitudes as function of plasma split parameter 
6. Arc current amplitude I a> 10 A; arc pulse fall time, T , 10 msec; K, 
5; total charge released, Q, 80 mC. 




Figure 8. — Transient amplitudes for time T in microseconds. 
Equivalent capacitance per cell, C /M 1 pF. 


Figure 6 shows the effect of the plasma split parameter 
5. <5=1 corresponds to C_ = l//?„=0 and 5 = 0, for 
which vj vanishes as it should, corresponds to 
C + = 1 /R+ =0. (See appendix A.) Amplitudes decrease 
gradually with decreasing 5 and then vary sharply when 5 
is small. This is because the equilibrium electron 
conductances are concentrated in a small region near V + 
where they are all assigned when 5 is large. Large values 
of 5 redistribute only the much weaker ion conductances 
and leave the conductance at V + constant. Only when 5 
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becomes less than about 0.08, which is the fraction of the 
array initially floating positive, do the large electron 
conductances at V + begin to shift to V _ and cause 
distributional effects to be seen in the transients. 

Figures 7 and 8 show results for two regions of T, the 
arc pulse fall-time constant. The linear variation with f () 
merely reflects that the solutions are explicitly 
proportional to I Q . In figure 7, and u ( j are seen to be 
nearly proportional to I () , with little dependence on T and 
C p . This is because T is much larger than circuit time 
constants (e.g., R + C + or C 0 R L ) so the circuit acts 
mainly like a series pair of resistors driven by / v (/ ) . In 
figure 8, however, where T is of the order of /? + C + , 
some dependence on Tis seen in n_. 

The values of v (i in figures 7 and 8 show consistently 
the relation v ( j/I a ~ 0.2 ohm. This means that, for the 
present array parameters, load transients equal in 
amplitude to only 8 percent of the terminal voltage are 
produced by arcs of 200-A amplitude. Such arcs need 
only last 0.5 msec to transfer 80 mC of charge from an 
array to the plasma. This amount of charge is expected to 
be stored on a 1000 by 1500 cell array operating at 508 V 
and floating electrically in LEO. 


Conclusions 

The model presented here is intended to provide a 
simple means for estimating the transient response of 
high-voltage solar arrays in LEO to negative bias arcs. 


An analytical formulation has been made possible for 
both the plasma interactions and the solar cell 
characteristics by linear approximations. The model is a 
first step toward more complex analytical and numerical 
analyses of solar array transient behavior. 

An especially needed addition to the model is an 
independently determined current /voltage characteristic 
for negative bias arcs. Introducing this to the circuit 
equations would permit determination of current and 
voltage transients free of arbitrary parameters. Without 
it, the model can calculate results only for given pulses. 
For the special conditions treated, pulses within the limits 
of the linear approximations of the model are not found 
to cause clearly unacceptable load transients. An 
estimation of the effects of stronger arcs awaits a 
nonlinear analysis. 

Even if arc current amplitudes are small enough to 
cause only negligible single transients, many small arcs 
might combine to form an essentially continuous 
discharge. This could clamp the negative end of an array 
at or near the arc voltage threshold, which is a lew 
hundred volts negative. The resulting enhanced electron 
collection due to snap-over at the positive end of a high 
voltage array could cause significant power losses (see 
Plasma Current Collection, p. 2). This effect needs 
further study. 


National Aeronautics and Space Administration 

Lewis Research Center 

Cleveland, Ohio 44135, September 19, 1984 
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Appendix A — Calculation of Plasma Interaction Components 


Calculation of Floating Equilibrium 

For a chain with as many as 1000 cells, the voltage 
distribution V along the chain may be approximated by a 
continuous function of position. If e(i) is the 
current/voltage characteristic of a cell and £ locates a 
point in the chain, then 


d 2 V _ ( dt \ / di \ 
dz 2 ~ V di ) \dz) 


for 0<z<N 


(Al) 


Let the point z Q divide the chain into ion and electron 
collection regions. Thus, V(z G ) is the plasma potential, 
here taken to be zero. For ion collection at z<z Q , L<0 
and di = -a/V dz in the ohmic approximation. For 
electron collection at z>z 0 , di= -a e V dz and K>0. As 
defined previously, a/ = zxl 0~ 8 and a^ = 2.4xl0' 6 
mho. If the slight variation of cell operating point along 
the chain is ignored, the cells have a common dynamic 
resistance, de/di = -p. Substituting these expressions 
into equation ( 1 ) yields 


d^V 
~dz 2 


-pajV= 0 


for z<z 0 \ 


(A2) 


-pa e V= 0 for z>z Q | 

With the condition that V be continuous at z a , the 
solution is 


V = 


zlX^sinh X, 


(Z-Z 0 ) 


V = A sinh \ e (z~z 0 ) 


for z<z 0 


for z>z 0 


(A3) 


where X e = (pa e ) \i= (pa/)l /2 , and A is arbitrary. 
Current balance is achieved by requiring that ion and 
electron currents be equal. Thus, 



(A4) 


(A5) 


Solving this gives 



(A 6 ) 


With 0 ^ = 2. 4xl0- 6 mho and a/ = 2xl0^ 8 mho, 
z 0 /N = 0.92. Thus, 8 percent of the chain floats above the 
plasma potential. For ram ions, with a r =10 “ 8 mho, 
Z 0 /N = 0.94. 

These results are unaffected by increasing the number 
of chains from 1 to M. With M chains, a factor of M 
would multiply both a, and a e in equation (A 6 ), leaving it 
unchanged. 


Plasma Interaction Distribution 

Let 6 be the fraction of an array whose plasma 
impedances Z p are added together in parallel at the 
positive terminal K+. Thus, 5 is the plasma split 
parameter (see p. 4). Let 7 = s 0 /N. Then the following 
lumped circuit elements obtain at the terminals in the 
mixed-case approximation: 


Capacitance: 

C+=&NMC p 
CJe (1 - o)NMC p 


Resistance: 

For 5 < 7 , 

=oNM e 

T- = NM\ (7 - b)a e + ( 1 - y)a. 

For 5 > 7 , 

T- =(1 - 5)NXtotj 

R __ 



which, using equation (A3), yields the condition 
cosh A/z a = cosh \ e (N-z 0 ) 
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Appendix B— Circuit Analysis 


The circuit of figure 3 obeys the following set of transient 
equations with initial conditions t? + (0) = n_(0) = 0: 


dv+ dv„ 

a l - sr+ a 2 v t -c- sf -fv^O 

(Bl) 

f dv _ , dv + . 

h \-jj~+b 2 V- ~ C ~cir'f V+ =ls ^ 


where 


ai = C f +C 0 c - C 0 

&i=C_+C 0 f=±- 



K 


, „ PqRl 
Po + Rl 


b 2 =- 


\ 1 

R~Z + W 0 


If the Laplace transform 


q(s) = ilv(l) = e~ sl v(()dt 
J 0 


is used, equations (Bl) form a pair of first-order 
algebraic equations in the variable 5. The solutions for 
q + and q __ are 


q + (s)=p{s) 


cs+f 
D ’ 


?..($)=/?(*) 


^15 + ^2 

^ D 


(B2) 


where 
p(s) = £/,(/) 

D = (!) is + b 2 ) (a \s + a 2 ) ~ ( cs +f) 2 


i s {t) = I 0 {e~ t/T ~e~ Kt/ T ) for K > 1 (B4) 

Then, doing the integral in equation (B3) gives 
v + (t)= (B\e-* K/T + B 2 e~ t/T + + B 4 e^ 2 ^ 

where 


( Xl + J r) ( X 2 +J r) 

b 2 = 

(x,+ 

B 3 = (K- 1) 


fi 4 = -(K- 


The expression for v is identical but with c* and / 
replaced by and 02 • From u + and e_ one can easily 
construct the load transient v+—U-> 

Let us consider a simple but useful limit to the previous 
solutions. Suppose that the solar array internal 
capacitance C a is so small that it can be neglected and 
that all the plasma impedance is at the positive terminal, 
i.e., 6= 1. Then, with C_ = \ /R _ =0and the limit C 0 — 0 
taken last, one finds 



The denominator D may be put in the form 


D = A{s-\ l )(s-\ 2 ) 


v,U) = I {9 R 



£> kt / 1 

(Kt t - T) 


(> t/T 


where A-b\U\~c 2 and the X’s are solutions of the 
equation D = 0. The transforms for q + and q may be 
inverted using well known theorems with the following 
results (ref. 9): 


( T t -T)(Kt ^ -T) _ 


036) 


t’-t - [1//I(^i ~ Xj)] i [(c\| +y)exp(X,T) - (tX 2 +./) 

x exp^r)] / v (/ - t)(/t (B3) 

and v is the same but with t* and /replaced by a \ and a 2y 
respectively. These are the solutions for arbitrary / s (/) 
provided / v (0) = 0. Let us take the special case 


v , (/) - v (/) = I 0 R 0 (e {/r e Ki/T)= - R () i s (t) 

where t , R k C , . This expression for the load transient 
(tq v ) is to be expected since, with C 0 = (), the arc 
current flows through a pure resistance R () . In actual 
practice, C 0 is quite small and equation (B6) gives a 
rough estimate of transient amplitudes. 
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cell characteristics, plasma interactions, and power system features are modeled 
approximately by a linear, lumped element transient circuit, and the time domain 
equations are solved. Exact numerical results for solar array common-mode and 
load voltage transients are calculated for typical conditions. Acceptable load 
transients are found for a range of arc current amplitudes and time constants. 
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